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ABSTRACT 

It  is  common  to  estimate  a  distribution  by  means  of  a  step  function. 
Such  estimates  can  be  made  continuous  by  connecting  the  points  of  the 
steps  with  straight  line  segments.   In  this  paper  the  best  estimator 
of  this  class  is  found  for  data  which  is  uniformly  distributed  using 
minimum  risk,  then  this  risk  is  compared  with  those  of  the  sample 
distribution  function  and  the  Pyke  estimator. 
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I  ..  INTRODUCTION 

In  reference  \_\"\   ,  Prof.  Read  studies  the  properties  of  some  estimators 

i 

for  a  continuous  distribution.  In  particular  ,  he  computes  the  risk  function 

(using  the  assumption  of  the  uniform  distribution)  of  the  estimator  that 

r  th 

connects  the  points  (  X  . — r  )  ,  "where  X  is  the  r   smallest  observation  . 
*  v  r'n+1  '   '  x  * 

with  straight  line  segments  .  This  estimator  is  shown  to  be  better  than 
the  sample  distribution  function  in  an  asymptotic  sense  •  It  is  called 
the  Pyke  estimator  . 

It  is  natural  to  ask  the  question  ,  "  "How  well  can  a  distribution  be 
estimated  by  a  function  that  connect  the  points  (X  ,C  )  with  straight  line 
segments?"  The  purpose  of  this  paper  is  to  find  the  sequence  \  C  \  that 
minimizes  the  integrated  risk  function  (using  the  uniform  distribution) 
and  to  compare  that  risk  function  with  those  of  the  sample  distribution 
function  and  the  Pyke  estimator. 


IIJJEEERMIHATIGH  OF  THE  OPTIMUM  COEFFICIEHT3 

Let  <C«,Cof   .   .  •  ,C  .{  be  an  increasing  sequence  with  C  j=  Q  and 
^  1  c.  n+ 1  j  U 

C  .  =  1  .A  continuous  function  estimator  for  F(x)  can  be  defined  by 

IWrl 


(2.1)  H(x)  -  C  -¥ 


x-X 


X  ,-X 
r+1  r 


(C   -C  ) 
r+1  r' 


for  ;  X  ^=x— X  .  r  -  0,..,n 
'     r      r-f1 


where  X.^-X-  .  .  .  £-X    are  the  order  statistics  of  a  random  sample 
from  an  absolutely  continuous  distribution  function  F(x).  It  is  convenient 
to  assume  that  the  population  sampled  is  bounded  and  contained  in  a 
finite  interval  T  a,bj  .Thus,  we  can  define  XQ  =  a  ,  X  «  -  b  and  the 
risk  function  (R)  by 


(2.2) 


b  r  -.2 

R  =  E  J    F(x)-H(x)l   dF(x) 


For  convenience  the  interval  £a,b3  can  be  taken  to  LO,lj  .The 

estimator  K(x)  is  defined  piecewise  according  to  v/hich  of  the  random 

intervals  (X  ,X  . )  contains  x  .  Assuming  the  data  come  from  a  uniform 
r  r«vi 

distribution  The  joint  p.d.f .  of  any  two  order  statistics  .  say  X  and 

r 

X  .  ,  is  as  easily  espressed  .  Using  X  =  u  and  X  1  =  v  the  joint 
p.d.f.  is 


(2.3)    fT  T   (u,v)  =  ■ ur"1(1-v)n"r"1   forO^r^n 

V  T4-1       (r-1)!(n-r-1)! 


and  0  ^  u  ^  T  ^  1 

(2.4)  fT  Y  (0,v)  -  n(1  -  v)n"1      for  r  =  0  and  u  =  0 

,.c,    -      /  4N  _   n-1         for  r  =  n  and  v  -   1 

(2.5)  fx  x   (u»1)   nu 

n'  rn-1 

Let  define 

(2.6)  A=C    .c 


r    r+1    r 


So  that 


<2'7>  Cr=?,A3 


rj 
j=-0 


and  for  simplicity  f_  T   (u,v)  =  f  (u,v) 

a  .a  4        r 
r'  r+1 

Then  rewriting  (2.2)  again 


(2.8)     R. 


2 
=  Z      1  J  J    *-C A. I  f  (u,v)dudvdx 

r=0     O^rU^X^rV^I  L         V  -  U    J    r 


Now,  the  problem  is  to  find  the  G .  values  that  give  the  minimum  risk. 
Using  classical  optimization  technique  with  Lagrangian  form  the  problem  is  to 


n 
Minimise  $  a  E  -  S  (  X  A  .  -  1 ) 

n 
Subject  to  21  A*  -  1 

3*0  J 


using  X  for  the  Lagrange  multiplier  .  oince 


I  n       f  (  {  r      r_1        z  -  u    n2 

C2.9)    \^T  )))     U-ZA A   f  (u,v)d 

r-0  0^u^v-*1  L    j=:0  J   v  -  u  ^J  r 


udvdx 


-MIA  -  1) 
0-0  J 


The  partial  derivative  of  \  with  respect  to  A  is 


a$ 


k-1     x  -  u    x  -  u 


\\ 


j*0   °    V-  -  U   J  V  -u 


f,  (u,v)dudvdx 


n 

r=k+1 


r-1        x  -  u' 
c-ZA  -A  

;p0   J        V  -  U 


f  (u,v)dudvdx  +• —  =  0 
r  2 


for  k  -  1  ,  .  .  .  ,  n  -  1i 

The  partial  derivative  of  (j>  with  respect  to  &    and  using  (2.4)  is 


a$ 


atf         r    r  x  x  < 


dvdx 


z  \  \  Ux  -  z:  a,  -a  —  )f  cu,v) 

r=t    J  J  J  >0     3         r  v  -  u      r 


dudvdx  +  —=0 


The  partial  derivative  of  <$  with  respect  to  A  and  using  (2.5)  is 


aft      r  c  r       n-1  x-uix-u  „  , 

aA     J  J  L      j-o°      ni-uJi-u 


dudx  +•  — -  0 
2n 


Let 


f     f 

g(u,v)  * ,         Ir(D  =  ^4  fr(u,v)dudv 

v  -  u  ^Y*1 

W  2  ({      2 

Ir(g)  =  o£u^  g(u,v)fr(u,v)dudv     ,     Ir(g^)  =  Qiu:4  g"(u,v)fr(u,v)dudv 


Using  these  notations  ,   then  the  system  of  derivatives  is 


a$         1  k-1  1  1  n        J 

(2.10)    *  /    xL(g)dx-:FA      /    L(g)dx-A      J     LCg  )dx  +  Z.        Ul  (1)dz 

a^O*  ;p0     J  0       *  S       ^  r=k+1   0J      r 

n         r-1  1  n  1  \ 

-     Z       T.  A     /    l(l)cbc-     ZA      /    l(g)dx+-=0 
r=k*1     j-0     J  0  r=k+1         0  2 


for  k  -  1    ,     6     .     •     ,  n-1 

All  terms  in  the  above  eguations  reduce  to  Beta  functions  .  The  relationship 
between  Beta  and  Gamma  function 


J   „«       „  1     r(m)PCn)   (m-1)!(n-1)i 

BCm,n)  =  J  v^d  -  v)n"1dv  =  — = 

0  \     (m+n)     (m+n-1 ) I 


is  used  heavily.  7/e  proceed  to  evaluate  these  terms  . 


Evaluation  of 


/  xl^gjdx  = 


k-1 , ,      vii-k-1 
z u       11-v) 

v  -  u 


nl 


dudvdx 


(k-1)l(n-k-1)i 


nl 


Ck-1>lCa-fc-1) 


.\fI&V"~ 


dudv 


2    )     v  -  u 


nl 


1     V    t       1 


k+2, 4         Nn-k-1 


(k-1)i(n-k-1)l       3        k    0 


J    v       (1  -  v)  dv 


f1     k-f2,«  xn-k-K  }     k+2,«  xn-k-U     \ 

)    v       (1  -  v)  dv  +  —       j    v       (1  -  v)  dv  I 


k*1  0 


k+2  0 


1  1 

—    r    k+2< 


kfl     0 


f  v^U  - 


xn-k-1,  *        }     k+2 

v)  dv  +  — 


k+2     0 


(   v^Cl-v)n"^1d 


.] 


(2.11) 


3k*4 


6(n+1)(n+2) 
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Evaluation  of 


J1  xl  tOta*    (    (    [zur-1(.1-vf-r-1 

i  r  J    J    J  (x-1)!(n-r-1)! 


duclvdx 


?: —  \  f  C  (   2     2.   r-1.      vii-r-1,    . 

=    1  \  \   (v  -u  )u       (1-v)  dudv 

(r-1 ) ! (n-r-1 ) !        2     lJ  J 


nl  1  (1        A         0  , 


(r-1)! (n-r-1)!        2   [r     0 


1         A 


(2.12) 


rf2     0 


r*1 

(n+1)(n+2) 


J1  vr+2(1-v)n-r'ldv   t 


Evaluation  of 


r  f  f  x"u   k-1,,   vn-k-1       n!.  ...  .  ,  A 

— — —  dudvdx 


/  I^dx^l     uk"1(1-v) 


0  J  J  J  v-u  (k-1)!(n-k-1)» 


(k-1 ) !  (n-k-1 ) !  [  2  J  J  v-u 


v)     dudv 


-JJuk(1-v)n 


cuav 


11 


n! 


1  r  1       A 


(k-1)! (n-k-1 )!  I     2 


l1    I 

«-  k      0 


k+1 , ,      » n-k-1 , 
v        (1-v)  dv 


1        J 


—  /    vk+l(i-v)n-k-1avl     -_    f 
k*1   0  J  3*1   0 


k-t-1,.      s  n-k-1, 
v       (1-v)  dv 


(2.15) 


1 


2(nf1) 


Evaluation  of 
1 


I 


Ik(g*)dx  = 


k-1/4      N n-k-1 


.v  -  u 


u~   '(1-v)J 


(k-1 )! (n-k-1 ) ! 


dudvd:c 


n! 


(k-1)! (n-k-1)! 


(1    f-r    3       3 

t5  JJ(v-u)2 


-jj 


2       2 

V  -  u 


(v  -  u  )' 


k/.  »n-k-1,    , 

u  (1  -  v)  dudv 


J  J(v  -  u)2 


k+1,.  xn-k-1,    , 

(1-v)  dudv 


n! 


1(1        1 


(k-1)! (n-k-1)!        3    [k 
1  1 


i 


k*1,«      vn-k-1 


—      ]     v        (1-v)  dv 


k+1     0 


/kt-1/4  *n-k-1, 

v        (1   -  v)  dv 


i 


(2.14)       = 


30*1) 


12 


Evaluation  of 

/'  Ir(1)ta  =  Jj{ur-1(1-vf-r-1.      D! 


dudvdx 


(n-r)l(n-r-l)! 
r+i/4   vii-r-1 


(n-r)l(n-r-l)!   r  0 


(      r+i,4   vii-r-i, 
—  J  v   (1-v)     dv 


1   1 

f   r+-1/.   vii-r-1, 
J  v   (1-v)     dv 

r+1  0 


(2.15) 


1 
n+1 


Evaluation  of 

1  1 

I    I  (g)dx  is  same  as  evaluation  of  (    L  (g)dx 


(2.16) 


0  0 

1  1 


Jlr(g)cbc  = 


2(n*1) 

The  use  of  these  evaluations  (2.11),  (2.12),  (2.13),  (2.14),  (2.15), 
(2.16)  in  (2.10),  then  the  derivative  with  respect  to  A,   is 

3k+4   1  k-1     1 


(2.17)     _£$_- JTA  -_d 

2AV   6(n+2)  2  r-0  J  3     k 


1    n  n    r-1     1   n       X  (n4-1 ) 


(m-2)  r-k+1        r-k+1   j«0  °   2  r*kf1 


13 


3k+4  3(k+2)(k+1)      n 

3C  -  2(C    -  C  x  4  3(n+1) 6  21   C 

n+2    k     k*1    k)  (n+2)       r=k+1  r 


-  3  (1  -  C   )  +  3  A  (ii4-1)  -  0 

or  after  evaluation  the  derivative  (2.10)  is 

[0    im  J£  v  5k  +6k+2 

^•15'        C.  -  C,  .  +  6  51   C   "=  3*(n+1)  4  3n  - 


k    k+1      -V  4  r 

r^k+1  n+2 

for;  k  -  1,  .  .  .,  n-1 

The  derivative  with  respect  to  A.  is 
(2.19)  9$ 


-_JJ(I        A)     ntl„T)»-'    dTto 


3AQ  v         v 


n  r  f  r       £J  x  -  u  x 

+  Z1       I  Ux  -  2.A.  -A-  )f  (u,v)dudvdx  +  ~-  0 

r=1  7J  '  j-0  J  v  -  u    r  2 

c*J  j|2(l-v)n-1   dvdx  -  n^j  [(rf)2  (1-v)n-1dvdx 

+     21  (xlp(1)  dx-  S.    C    [l(1)dx-  lAr  fl(g)dz44=  ° 
rM  '  r=-1        '  r»1       ' 


This  can  be  reduced  to 

n  2 

(2.2Q)      -U6I  C.s  3>(n+1)  +  3n  -  -^-       f  or  ;  k  =  0 
1     k=1  *  n+2 


The  derivative  with  respect  to  A  ig 

n 


(2.21)     Q$    f  fr    ^ 


dudx  +  ~  -  0 


14 


can  be  reduced  to 


ajj        2 

= A  -  3A(n+1)  -  0 

dA        n+-2        n 


or 
(2.22) 

so  in  general 

(2.23) 


Cn=  3Mn+1)  4 


n+2 


for;  k=-n  . 


Cjc  -  cwi  +6     A     C,-  3Mn+1 )  4  3n  - 


k«*1  -r~«     r 


3k  +6k42 


n42 


for  ;     k  =  0,   .   .   .   ,  n 


Using  (2.23)    ,   then  system  of  eguations  in  matrix  form  is 


5  6  6  6 
15  6  6 
0     15     6 


.  6  0 
.  6  0 
.     6     0 


0     0     ...     1     5     6     0 


a    0     .     .     ♦     0     1     5    0 


0     0     .     .     #     0     0     1  -1 


where  6  s  3A(n+1)  4  3n 


r«, 


n-1 


n 


n+1 


6  -  2/(n  4.  2) 
e  -  11/(n  4  2) 
6  -  26/(n  +  2) 


e  -  (3n  -6n42)/(n+.2) 
6  -  (3n2-1)/(n+2) 
6  -  (3n246n42)/(n42) 


15 


After  successive  row  subractions  the  matrix  is 


4     10     0 
14     10 

0     14     1 


0  0.. 
0  0.. 
0     0.. 


..00 
..00 

..00 


14  10 
0  14  1 
0     0     1-1 


3 


9/(n+2) 

15/(n+2) 
21/(n+2) 


i  Cn-1   ! 


C 
» 

LCn+1 


(6n-3)/(n+2) 
(6n+3)/(n+2) 


j      |_ 


|6-(5n  4.&H.2)/(n+2) 


To  solve  this  matrix  equation  i3  considered  the  linear  second-orde: 


difference  equation  such  that 


(2.24) 


S  ♦  4(W+  <W  - 


Gk  +  9 


n  +  2 


with  initial  conditions  Crt  -  0  and  C  .  =  1 

0        n+1 

The  corresponding  homogeneous  equation  of  (2.24)  i3 

(2.25)  \  -  4Cbfi  4  0k+2  =  0 
has  the  auxiliary  equation 

2 

m  +  4m  *  1  a  o 

•Che   roots  of  this  equation  can  be  calculated  using  the  quadratic  formul: 
These  roots  are  A  =  -2  + 15  and  B  -  -2  -  \T5"1 
so  the  general  sotion  of  (2.25)  is  given  by 

(2.26)  C^  =  K1Ak  4  lyi*  +  C* 


16 


* 
Trial  solution  of  C.  is 

k 


i 


Using     (2.24) 


£*D0  +  D1k 


Drt  =     „/1   „v    *  D4-  — Kr         then 

0         2(n-v2)  1     n+2 


*  _    1   4  2k 

°k" 


2(n+2) 


So  the  general  solution  is 


(2.27)         Cfc  ~-  K,Ak  +  K/  ,  $*$- 

using  initial  conditions  to  determine  the  values  of  E.  and  Kp 


CQ  «  0  in  (2.27)  implies 
(2.28)  K.  f  K0  -  - 


1  "*  "2     2(114-2) 


C  «  a  1   in  (2.27)  implies 
n+1 


.n+1   _  _n-*-1  __   1 


(2.29)  V      ^H/     -aisar 

Simultaneous  solution  of  (2.28)  and  (2.29)  are 


1  t  B^1  It  A1*1 


K  -  __ — — _ __  aad  r   -  _ 


1  "  2(n+2)  (An+1  -  Bn+1>       2    2(n+2)(An*1  -  Bn"  ) 


so 


(„„,        ^-LL^r-^i-^-^) 


(nr2)  I  2(Anv1  _  Bn+1} 


for;  k  -  0,  .  .  .  ,  n+1 


17 


usiiig     \  =  C         -  Ck         then 


(2.51)       A  1        (  (1   -  A)(An"k4  Ak)  -  (1  -  B)(Bn'k+  Bk) 

k         (n+2)   [  2(AD*1  -  B^1) 

for;       k  =  0,   .    .   .   .   ,n 


18 


Ill, CALCULATION  OF  RISK 
Prom  equation  (2.8)  the  risk  equation  is 
(3.1)  »«EJ]j[(x-Cr)-  k?Aj|  fr(u,v)dudvdx 


n  r  •  1     ?  J  ~     1 

21  \    J     xV(1)dx  -  2C     J     xl  (1)  dx  4-  CM  I  (1)dx 
r^OO  rOr  rOr 


-  2  A     J     xl   (g)dx  +  2C  A     J     I   (g)dx  ^A^     /     I   (g^)dx 
r  0  rr0  0 


Using, 

f1  -A[in«.-i  istra* , 

'  rN  (n4-3)(n4-2)(n4-1 


j     xlr(1)   dx  = 


(n+1)(n+2) 


0       r  n+1 


0         r  6(n+1)(n+2) 


/   V^ta      =    ^TnTT) 


r  =  0,    .    .   .,n 


1 

r+1  r  =  0,    .    .    .,n 


f    i   0)dx      =  —. '  r  =  0,    .    .    .,n 


(    xl  (g)dx    ^    3r  +  4     ..  r  =  0,   .   .   .,n 


r  =  0,    .    .    .,  n 


1       '  O  1 

f    I  (e-  )dv    =    .  r  =  0,    .    .    ,,n 
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So 


3         3ln«-2)         Tn*-2)  £-j       r         n4-2    ^-« 


n  ,       n 

C 


«     n  n 

5     r=0  r=0 


*L      2 
Calculation  of     x    A 
r5)     r 


Using  (2.31) 


±4     =    -i-2  Z  5   1   -         n+1      1   SI    [(^)Un-r  +  Ar)  -   (1-B)(B^4  A 

fe)  r  (^-2)^0    L  2(An+1   -  Bn+1)L  JJ 


and  using 


1  -  An+1 


=  2 


T         (A*-r  +  Ar)  =     2T    i 

—  ^~  1    -  A 

r=0  rrO  1       A 


Z  (Bn~r  t  Br)  =  2Z 


n 

^  =2 


1  -B11*1 


r=0  r-0  1  -  B 

Jimiliarly 


r=0  r=0  1  -  A2 


20 


1        B2(m"1) 
Z     (D2(n-r)   +   B21   )     =  2   Z    B2*     =     2   — 1^ 

r=0  r=0  1  -  B 


and 


1  1    -  A 


\<m  A"2         1  -  B2         6(14B) 


1  _  1         ^      1    -  B 


1   _  B  "2       r«  A2         6(1   4  A) 


1   -  A     _        1   -  B 


1  4  A  1   +  B 


B' 


[L    (1   -A2(n+1))     »-±-   (Bn+1   -An+1) 


1  +  A  1  4  B 

then 

,,0,        <§-  A  2       n*3  xm-1  (1  -A)2An  4  (1   -  B)2  B*  -  12 

(5.2)       2_AT.=^ — ~2     +  2 


rs 


-0  (n+2r         (n+2r  2(An4l  _  ^a+1j  2 


-f? 


(An+1   +Bn*-1)  +  2 


2(n+2)2(An+1-Bn+1) 


Calculation  of   7      c     1   c 
r=0 


Using     (2.30) 

1 


C„  ,   C_  ^    *  ^     \     (r2  4  2r  +  |)  4  (r  ^J^C***)^   -r-|] 

r-0        V- 


2-     »1     '         (n+2)2 
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+  (r +-|)  [(o*2)  Or-r--ij 

1  U-1-*  -  a-  -  *»'-  ♦  rf  [  a-  -  e*  .  rf« *  b~'J 


4 


4  (An4l_Bn4l}2 


here 

T     (r2  4  2r  +-2)     -     n(n+1)(4n+U)  4  9(n*1)_ 
r=0  4  12 


and 


^  _±_  [in4"1"r  -  Ar  -   ^UT  4-   Br][  An"r  -  A1*1   -   rf"  ♦  B^1] 


r=0  4  (AnM  _  Bn4-1    )2 


^t2*"3  *  B2n+5  '   (6n4-4)(An+2  4  Bn+2)  -   (A21141   4-  B^1  ) 


.(fiw.7)(An  +  Bn)*8(,rt-1U 

Then 

(3-5)     icr+,Or-    .    n(n„)(4n+14)  ,  9(nH)  4jT     ,|+  2£  £     1 

r^O  r+1   r  12(n+2)^  n+2  L         *        r«1       r    r*1     rJ 

[a^V**5  -(6n+4) ( A^+B^J-CA20*'1^20*  V(6n+7) (An.Bn)+  8<in 

24(n+2)2(An+1   -B*"1)2 
Using  (3.2)   and   (3.3)      then 

2n2+.3n-1         \T?  (A™"1    ^  ^  ) 
4)      (nv1)  R^- 


12(n*2)2      "  12(n+2f  (An+1-?nT 

Li  j*1 ,2  -[(2°*3)(A\^)  4  40  (n+1) 


24(n+2)2(Axl-r,-B1"-') 


42n(An+2  4  Bn+2)   4  8(a+1)    (A^1   4  BnM) 
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1 


ir»COHPAHISC!fS  AND  simliary 

Prom  (3.17)  page  14  of  Prof .  bead's  paper 

(n+1)  R(C)  =  J  (n-1)x(1-x)dx  +  2(n4-1 )  J  \    (1-x)h(x)  +  xh(1-x)(  dx 
0  0  L  J 

where 

n 


r  x  /  x  -  u 

n+t 


C  = and       h(x)  =  /  !  j  udu 

0   M  -  u  i 


2  *   n~1  /-1 

(ml)'TlCC)  = +  4(n+1)    (1-x)b(x)dx 

6  0 


and  evaluation  of 


1  1  x/x-u  \ 

[  (1-x)h(x)dx  =  /  (1-x)  /[— — lududx 

0  0  0  \1-u  / 

=    I  I  Cl-x)(—  \w 
O^x^-I      v  1-u  ' 


let 


X-U 

y  =  _  Q^y^ 

1-U 

then 

dx 
X  «  XL  +  (1  -  u)y  ,  1  -  x  =  (1  -  u)(1  -  y)  and  dy  =  


1  -  u 


so 


n-1         1 


(4.1)  (n+1)R(C)  -  + 


6(n+.1)   3(n+2)(n4l) 
where 

R(C)  is  Pyke  risk 
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Risk  of  sample  cumulative  distribution  function  R(SDF)  is 


1        t  1 

R(SDF)  -—      I     x(1-x)dx  m  — 
n     0  6n 


then 

n+1 
(4.2)  (n+1)R(SDF)  =  

6n 

Optimum  risk  R(OPT)  is  derived  in  this  paper  is 

2n2*3n-1  f5{k*+KiP+y «) 

14.3)  Cn4l)RC0PT)  -  . ir~     ~   2~"n7r~a7l 

12(n+2)2      12(n42j^(An^1-rf^1) 

(2Bf 3) (AW)+8(-m-1 )  (An+1.»Bn*1  )*2n(AXM'2^4'2)440(x>4.1 ) 

24(n+2)2CAn+1-Bn+-1)2 

Lot  this  eguation  last  two  terms  may  be  omitted. Since  they  are  close  to 

zero  for  n  —  6 

so 


(n*1)R(0PT)  = 


2n2+3n-1 


12(n+2)2 


Cumputer  output  is  shovm  on  page  25  for  n  =  1  ,  •  •  .  ,  60  .  Equations 
(4.1)  ,  (4.2)  and  (4.3)  are  used  . 


24 


N 


(N+1)R(C) 


(N+1)R(SCF) 


(N+1)R(0PT) 


LAST  TWO  TERMS 

OF 
(N+l)R(OPT) 


1 

0.05555555 

2 

0.08333331 

3 

0.09999996 

4 

0.11111100 

5 

0. 11904750 

6 

0.12499990 

7 

0.  12962950 

8 

0.13333320 

9 

0.13636360 

10 

0.13888870 

11 

0.14102550 

12 

0.14285700 

13 

0. 14444440 

14 

0.14583320 

15 

0. 14705870 

16 

0.14814800 

17 

0.14912270 

18 

0.14999980 

19 

0.15079360 

20 

0.15151500 

21 

0. 15217380 

22 

0.15277770 

23 

0. 15333320 

24 

0.15384600 

25 

0.15432090 

26 

0.15476170 

27 

0.15517230 

28 

0.15555540 

29 

0. 15591390 

30 

0. 15624980 

31 

0.156  56550 

32 

0.15686260 

33 

0.15714280 

34 

0.15740730 

35 

0.15765750 

36 

0.15789460 

37 

0. 15811950 

38 

0.15833320 

39 

0.15853650 

40 

0.15873000 

41 

0.15891460 

42 

0.15909080 

43 

0. 15925920 

44 

0. 15942010 

45 

0. 1595743C 

46 

0. 15972210 

47 

0.15986380 

48 

0.15999990 

49 

0. 16013060 

50 

0.16025630 

51 

0.16037730 

52 

0.16049370 

53 

0.16060590 

54 

0. 16071420 

55 

0.16081860 

56 

0.16091940 

57 

0.16101680 

58 

0.16111090 

59 

0.16120200 

60 

0.16129020 

0.33333330 

0.25000000 

0.22222220 

0.20833330 

0.19999990 

0. 19444440 

0.19047610 

0.18750000 

0.18518510 

0.18333320 

0.18181810 

0.18055550 

0.17948710 

0.17857140 

0.17777770 

0.177C8330 

0.17647050 

0.17592590 

0.17543850 

0.17499990 

0.17460310 

0.17424230 

0.17391300 

0.17361110 

0.17333320 

0.17307680 

0.17283940 

0.17261900 

0.17241370 

0.17222210 

0.17204290 

0.17187500 

0.17171710 

0.17156850 

0.17142850 

0.17129620 

0.17117110 

0.17105250 

0. 17094010 

0.17083320 

0.17073160 

0.17C63480 

0.17054250 

0.17045450 

0.17037030 

0.17028980 

0.17021270 

0.17013880 

0.17006800 

0.16999990 

0.16993460 

0.16987170 

0.16981130 

0.16975300 

0.16969690 

0.16964280 

0.16959050 

0.16954020 

0.16949140 

0.16944440 


0.05353007 

0.07628661 

0.09249049 

0.10353520 

0.11179110 

0. 11814010 

0.12318110 

0.12727660 

0.13066930 

0. 13352540 

0.13596240 

0.13806620 

0.13990060 

0. 14151430 

0. 14294460 

0.14422110 

0.14536740 

0.14640240 

0. 14734160 

0.14819760 

0.14898100 

0. 14970070 

0. 15036420 

0.15097770 

0. 15154670 

0. 15207600 

0. 15256950 

0. 15303060 

0.15346260 

0.15386800 

0. 15424930 

0.15460860 

0. 15494760 

0.15526810 

0.15557160 

0.15585930 

0.15613250 

0. 15639220 

0. 15663940 

0.15687500 

0.15709970 

0. 15731440 

0. 15751970 

0. 15771620 

0.15790440 

0. 15808480 

0. 15825800 

0.15842430 

0.15858420 

0.15873790 

0. 15888600 

0. 15902860 

0.15916600 

0. 15929870 

0.15942670 

0.15955030 

0. 15966980 

0. 15978530 

0.15989710 

0.16000530 


-0.01649304 

-0.00857831 

-0.00582393 

-0.00399830 

-0.00294773 

-0.00225484 

-0.00178204 

-0.00144335 

-0.00119288 

-0.00100234 

-0.00085407 

-0.00073642 

-0.00064150 

-0.00056382 

-0.00049944 

-0.00044549 

-0.00039983 

-0.00036084 

-0.00032730 

-0.00029822 

-0.00027285 

-0.00025059 

-0.00023094 

-0.00021352 

-0.00019799 

-0.00018410 

-0.00017163 

-0.00016038 

-0.00015020 

-0.00014095 

-0.00013254 

-0.00012486 

-0.00011783 

-0.00011137 

-0.00010543 

-0.00009996 

-0.00009490 

-0.00009021 

-0.00008586 

-0.00008182 

-0.00007806 

-0.00007455 

-0.00007128 

-0.00006821 

-0.00006534 

-0.00006265 

-0.00006012 

-0.00005773 

-0.00005549 

-0.00005338 

-0.00005138 

-0.00004950 

-0.00004771 

-0.00004603 

-0.00004443 

-0.00004291 

-0.00004146 

-0.00004009 

-0.00003879 

-0.00003755 
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After  computer  calculations  for  n  =  1  ,  .  •  •  ,  60  (4.1)  ,  (4.2)  and 

...    I 

(4.3)  is  shown  in  (.4.4)   and  (4.5) 


(m-1  )R 


.35 


.30 


.25 


.20 

.16 

.15 


.(n+1)R(3DF) 


10 


.05 


10 


20 


30 

(4.4) 
26 


40 


50 


» 

60 


E 


.050 


.025 


.020 


.015 


.010 


.005 


10 


20 


30 


40 


50 


60 


(4.5) 
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